Long range magnetic ordering in Na2lr03 
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We report a combined experimental and theoretical investigation of the magnetic structure of the 
honeycomb lattice magnet Na2lr03, a strong candidate for a reahzation of a gapless spin-liquid. 
Using resonant x-ray magnetic scattering at the Ir La-edge, we find 3D long range antiferromagnetic 
order below Tjv = 13.3 K. From the azimuthal dependence of the magnetic Bragg peak, the ordered 
moment is determined to be predominantly along the a-axis. Combining the experimental data 
with first principles calculations, we propose that the most likely spin structure is a novel "zig-zag" 
structure. 



PACS numbers: 75.25.-j, 75.30.Et, 75.30.-m 

Recently, Ir-based transition-metal oxides have at- 
tracted significant attention [lUlOj . owing to a combina- 
tion of the strong rclativistic spin-orbital coupling (SOC) 
of the Ir bd electrons and the non-trivial topology of 
underlying lattice in many of the structural families. 
These latter include the pyrochlore A2lr207, the hyperk- 
agome structures of A4lr308 and the honeycomb lattices 
of A2lr03, where A is an alkaline metal or lanthanide. 
These offer the possibility for exotic electronic behavior. 
Indeed, the strong SOC in the iridates has been shown 
to lead to novel Mott-insulating states even in the pres- 
ence of only weak correlations . The non-trivial under- 
lying lattice geometries enhance the significance of this 
coupling because the spin and orbital components of the 
wavefunction then become entangled, and this can lead 
to various kinds of strongly anisotropic magnetic interac- 
tions, depending on the topology of the underlying lattice 
[2]- Thus these materials provide an array of interesting 
physics, with exotic states such as topological insulators 
[31 E] and spin-liquid states ISHTU] proposed as possible 
ground states. 

Among the iridates, A2lr03 is particularly interesting. 
Here ions reside in oxygen octahedra [TTJ [T^] , and 
are suggested to be magnetic with an effective moment, 
Jg^^ = 1/2 gmillll]. This combination of Jeff = 1/2 
and the underlying honeycomb lattice makes A2lr03 a 
promising candidate for the long-sought realization of the 
exactly solvable Kitaev model |T3j. This model describes 
spin S=l/2 magnets with extremely anisotropic exchange 
on a honeycomb lattice. Its ground state is a gapless 
spin liquid, which has attracted much attention for its 
potential in quantum computing applications due to its 
insensitivity to local perturbations. 

Unfortunately, experimental studies of A2lr03 have 
been rather limited to date. For Na2lr03, the only ex- 



perimental information is the very recent work of Singh 
et al. who carried out susceptibility measurements and 
suggested that Na2lr03 is antifcrromagnetically(AF) or- 
dered at low temperature [12], apparently deviating from 
Kitaev model predictions . This raises the important 
question of why does the Kitaev model break down in this 
system? To address this issue, an essential first step is 
to understand the nature of the low temperature ordered 
state, including the magnetic structure, the ordered mo- 
ment direction and dimensionality of the ordering. 

Here we report x-ray resonant magnetic scattering 
studies of Na2lr03, performed at the Ir L3 edge. Our 
results explicitly show that Na2lr03 has a long-range an- 
tiferromagnetically ordered ground state below Tjv=13.3 
K and that the ordering is three dimensional(3D). From 
the azimuthal dependence of the magnetic peaks, the or- 
dered magnetic moment is determined to be mainly along 
the crystallographic a direction. Based on our observa- 
tions, two magnetic ordering structures (see Fig. 3) are 
found to be possible candidates: "zig-zag" and "stripy" . 
Combining the experimental data with a set of con- 
strained first principles calculations, we propose the zig- 
zag phase as the most likely ground state. 

Na2lr03 crystals were synthesized by a solid-state syn- 
thesis method. The synthesis and bulk property char- 
acterization are reported in Ref.[T2]. The x-ray reso- 
nant magnetic scattering experiments were carried out 
at beamline X22C at the National Synchrotron Light 
Source, Brookhaven National Laboratory. The single 
crystals used in our experiment were plate-like with typ- 
ical dimension of ^1x1x0.1 mm"^. A Ge(lll) analyzer 
was used to improve the resolution and to suppress the 
inelastic background due to fluorescence. Except for the 
temperature dependent work, all the data were collected 
at 1.58 K. 



2 



0.8 



0.6 



0.4 



0.2 



absorption 

• magnetic peak (0 1 11) 
— struoturai peal^ (0 2 1 2) 




11.145 11.175 11.205 11.235 

Photon Energy(keV) 



£2 4500 



1500 



T 1 r 

% 




-0.02 0.02 0.9 0.95 1 1.05 1.1 10.98 11 11.02 

H(r.l.u.) K(r.l.u.) L(r.l.u.) 

FIG. 1: (Color online). Top panel: Incident photon energy 
dependence of the (0 1 11) (circles), and the (0 2 12) (dia- 
monds) Bragg peaks on tuning through the Ir L3 edge. The 
solid line is the absoption spectra of the L3 edge as measured 
by monitoring the fluorescence signal. Bottom panels: H, K 
and L scans across the (0 1 11) peak. 



At low temperatures, well-defined peaks were observed 
at (0 ±1 L) Bragg points with L—7, 9 and 11 with the in- 
cident x-ray energy tuned to the Ir L3 edge. For Na2lr03 
which has a C2/c structure [12], these are structurally for- 
bidden Bragg points. In Fig. 1, we show the energy de- 
pendence of one of these peaks, the (0 1 11), near the Ir 
L3 edge, together with that of the structurally allowed 
Bragg peak (0 2 12). In contrast to the (0 2 12), the (0 1 
11) peak shows a large enhancement on tuning through 
the Ir L3 edge, as determined by the fluorescence sig- 
nal (also plotted). In principle, this resonant behavior 
could be due to either orbital ordering or magnetic or- 
dering. To distinguish the two is non-trivial, especially 
when there is strong SOC. Based on the fact that the po- 
sition and shape of this resonance is similar to that previ- 
ously observed in Sr2lr04 [IT and that the temperature 
dependence of the ordering is consistent with the sus- 
ceptibility measurements [T2] . we conclude that the most 
likely explanation for this peak is that it is magnetic in 
origin. That is the (0 1 L=odd) reflections demonstrate 
the presence of long range antiferromagnetic order. 

Scans through the (0 111) peak along the three high- 
symmetry directions are shown in the bottom panels of 
Fig. 1. The peaks are at, or close to, the resolution limit 
in all three directions, with the resolution determined 
by similar scans performed through the nearby (0 2 12) 



allowed Bragg peak. This indicates that the magnetic 
ordering is 3D in nature and is long ranged. From the 
half-width-at-half-maximum (HWHM) of these scans, we 
place lower limits on the respective magnetic correlation 
lengths, defined as {jjy^jjjf), to be about ^mag > 6OOA, 
120 A and 500 A along the crystallographic o, b and c 
directions, respectively. We emphasize that these values 
largely reflect the instrumental resolution and that the 
resolution along the K direction is the broadest. 

In order to determine the AF ordered moment direc- 
tion, the azimuthal dependence of the (0 1 11) peak was 
measured by rotating the sample about the (0 1 11) di- 
rection. In our experimental setup, the incident beam 
is a polarized, that is the electric field is perpendicular 
to the scattering plane. In this geometry, the dipole res- 
onant magnetic x-ray scattering cross-section varies as 
kf ■ M, where kf is the direction of the scattered beam 
and M is the magnetic moment direction [T5]. In Fig. 2, 
the intensity of the (0 1 11) magnetic peak is plotted as 
a function of a, where a is the angle between kf and 
the a-axis. The a dependence of the intensity follows 
cos^ a (dashed line) reasonably well, suggesting that the 
magnetic moment is predominantly along the a-direction. 
The slight deviations from this curve may indicate that 
perhaps the magnetic moment also has some small com- 
ponent in other directions. 
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FIG. 2: The azimuthal dependence of the (0 1 11) magnetic 
peak intensity, a is the angle between scattered beam and 
the crystallographic a direction. 

We also carried out a search for other peaks at various 
positions, and along high symmetry directions, in recip- 
rocal space. In particular, scans across (0 1 L) with even 
values of L showed no sign of magnetic peaks, indicating 
that the two Ir sublattice planes, which are separated by 
half a unit cell in the c direction, are antiferromagneti- 
cally coupled. Further, no peaks were found between (0 1 
11) and (0 11), or between (0 1 11) and (1 1 11). These 
latter observations, together with the determination of 
the in-plane ordering wavevector to be (0 1), allow us to 



3 



conclude that the magnetic unit cell size is the same as 
the structural unit cell. Unfortunately, scans across (0 3 
L) and (1 L) with L = 11 and 13 were contaminated 
by nearby strong non-magnetic powder peaks, which pre- 
vent us from solving the magnetic structure completely. 
Nevertheless, the observed peaks do provide strong con- 
straints on the possible magnetic structure of Na2lr03, 
as we shall now discuss. 
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TABLE I: The total energy Etot per Ir, the bandgap Egap and 
the absolute value of the spin | (S) \ and orbital j (L) \ moment 
of four coUinear magnetic configurations calculated within the 
LDA-I-U-I-SO approximation with the moments constrained 
along each of the three crystallographic axes (a, b and c). 
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FIG. 3: Crystal structure and possible antiferromagnetic or- 
dering patterns of the Ir sublattice. (a) crystal structure: the 
arrows at the left-bottom corner show the in-plane directions 
of the crystal; the solid-line box indicates the in-plane struc- 
tural unit cell; different Ir-Ir bonds are labeled with L(long), 
M(medium) and S(short) with the bond lengths to be 3.30A, 
3.23A and 2.86A respectively, (b), (c) and (d): 3 possible 
magnetic structures. In each case the magnetic unit cell is 
the same as the structural unit cell shown in (a). The shaded 
boxes highlight the stripy and zig-zag chain elements accord- 
ingly. 

With the determination that the magnetic unit cell 
is the same as the crystallographic unit cell and that 
the ordered moment points along the a-axis, the number 
of possible AF configurations in Na2lr03 are reduced to 
three, shown in Fig. 3. The Neel and so-called "stripy" 
states were discussed by Chaloupka et al.^U\ and were 
found by adding an isotropic Heisenberg term to the Ki- 
taev Hamiltonian, while the zig-zag pattern has not been 
considered in Na2lr03 before. The non-zero (0 1 L) peaks 
rule out the Neel state, leaving just the stripy and zig-zag 
patterns as possibilities. 

To further resolve the magnetic ground state, we per- 
formed spin density functional theory calculations 17J of 
four coUinear magnetic configurations (FM, Neel, stripy 
and zig-zag) in which the moments are constrained along 
each of the three crystallographic axes, using the em- 
pirical monoclinic C2/c structure [T^]. The results are 



summarized in table |T] For all the moment directions, 
the zig-zag configuration is found to have lower total en- 
ergy than the stripy configuration. This, in combination 
with the x-ray analysis, leads us to propose the zig-zag 
configuration as the most likely ground state of Na2lr03. 

The calculation results also show that for the in-plane 
directions of the moments (a and b) , the FM and zig-zag 
configurations are nearly equally low in energy, whereas 
the stripy and Neel configurations are nearly equally high 
in energy. The other electronic properties, the band gap, 
the spin and orbital moments, display the same trend. 
We observe that the FM and zig-zag configurations have 
ferromagnetically aligned moment in the shortest Ir-Ir 
bond, while the stripy and Neel have the same bond an- 
tiferromagnetically aligned (see Fig. 3). We note that 
the constrained moment directions investigated in this 
study are not meant to represent the moment direction 
of the theoretical ground state, which in general could be 
away from the crystallographic axes and non-coUinear. 

We now turn to the temperature dependence of the 
magnetic ordering, which was monitored by H scans 
across (0 1 11) peak. The magnetic peak, whose width 
is almost a constant as function of temperature, was fit 
to a Lorenztian-squared function. In Fig. 4, we plot 
the temperature dependence of the integrated intensity 
obtained from fitting. From the onset of the magnetic 
peak intensity, Tjv is determined to be 13.3(±0.1) K. We 
note that the order parameter could not be fitted to a 
simple power law, even for the region near T/y. There 
is a small linear region close to the transition which ex- 
trapolates to a lower T^r = 12.7 K(see Fig. 4 insert). 
Above this temperature, a long tail persists up to 13.3 
K. One possible explanation for this is a small rounding 
of Tjv arising from sample imperfections, though other 
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explanations are also possible in systems with complex 
magnetic interactions. For example, similar temperature 
dependence was observed in the 4/ heavy femion sys- 
tem CeAs|16j. There the gradual onset of long range AF 
ordering was attributed to the suppression of the devel- 
opment of long range magnetic ordering due to strong 
short range spin fluctuations. In Na2lr03, the presence 
of short range correlations have also been suggested from 
magnetic and heat capacity measurements [T^ . 
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FIG. 4: Temperature dependence of the integrated intensity 
of the (0 1 11) peak. Insert: Expanded view of the tempera- 
ture region near Tjv. 

We next discuss the significance of these results in the 
context of the search for an experimental realization of 
the Kitaev model. There is a general consensus that sys- 
tems with effective spin 1/2 on a honeycomb lattice with 
strong SOC, as are present in the Na2lr03, are good 
places to look. For example, Chaloupka et aL[Tn] showed 
that for such a system the underlying Hamiltonian could 
be described as a sum of an isotropic Heisenberg term 
and the Kitaev Hamiltonian. They found a range of pos- 
sible ground states, from a Neel state to the quantum 
spin liquid. Our data rule out the Neel state, implying 
that the Kitaev term is significant in Na2lr03 and that it 
is not governed by a simple Heisenberg Hamiltonian. A 
second point to emphasize is the observed magnetic or- 
dering is three dimensional, that is the correlation length 
is as long perpendicular to the planes as within them. 
The Kitaev model describes a two-dimensional honey- 
comb lattice and its ground state is expected to be robust 
to small perturbations. The fact that we are not seeing 
the spin liquid state suggests that these perturbations, 
which might be out-of-plane interactions, for example, or 
finite Heisenberg exchange, are relatively strong in this 
material. 

In conclusion, three-dimensional long-range antiferro- 
magnetic ordering is found in Na2lr03 by resonant x-ray 
scattering at Ir L3 edge. The azimuthal dependence of 
the magnetic scattering intensity indicates that the or- 



dered moment is predominantly along the a-direction. 
Two possible ordering patterns, namely stripy and zig- 
zag, are consistent with the experimental data. Com- 
bining the magnetic x-ray analysis with a set of con- 
strained first principles calculations, we propose the "zig- 
zag" phase as the most likely ground state. Given that 
no real materials have been found to show pure Kitaev 
physics and that Na2lr03 was promising from many as- 
pects, understanding our observations from a theoretical 
perspective could shed light on why the Kitaev model 
breaks down in Na2lr03, and thus pave the way towards 
the final realization of spin-liquid states in real materials. 
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